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Abstract 

We present a calculation of full one-loop radiative corrections, including the constant 
term, to the asymmetry parameter of polarised neutron beta decay. This gives the radia- 
tive correction to the axial coupling constant gA extracted from the beta asymmetry so 
that it ties to gA that appears in neutron decay lifetime in a consistent renormalisation 
scheme. We find that the ratio of axial-vector to vector couplings determined from the 
beta asymmetry, after taking account of the outer radiative correction, is related to the 
bare value as G A /G V = 1.0012GVG^. 
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Beta decay asymmetry of the polarised neutron has been used to determine the axial 
vector coupling constant gA of the nucleon. One-loop radiative corrections to the asymmetry 
parameter has been calculated by several authors CQ-jH, but their calculations do not include 
the constant term, or so-called inner corrections, which require a special care in the treatment 
of the UV divergence of the radiative correction. This makes the identification of gA extracted 
from p decay asymmetry with that which appears in the nucleon beta decay rate ambiguous 1 . 
The radiative correction to beta decay is UV divergent and it is rendered finite only with 
the use of Weinberg-Salam's theory of electroweak interaction. The complication arises from 
the fact that one must deal with quarks in electroweak theory, and one must continue to 
calculations with hadrons at low energies 0. 

In this paper we calculate one-loop radiative corrections including the inner correction: 
hereby the coupling constant qa that appears in the (3 asymmetry parameter is unambiguously 
tied to that in the decay rate in a consistent renormalisation scheme. The key point of the 
calculation is the clarification of the universal and non-universal UV divergent parts by using 
the current algebra technique and the proof that the same combinations of the renormalisation 
factors appear in the beta asymmetry as in the beta decay rate. The separation of the 
UV divergences into universal and non-universal parts was done first by Abers et al. jS] 
for the Fermi transition of nuclear beta decay, and then used by Sirlin |5] to develop a 
practical scheme of one-loop radiative corrections for the + — > + transition. The scheme 
was extended to the Gamow- Teller transition by our recent publication (Paper I) 8 . The 
present work is an application of the formalism developed in Paper I. We are content with 
the outline of the calculations in this paper, since the bulk of lengthy calculations are parallel 
to the ones presented in Paper I. We refer the readers who are interested in technical details 
of calculations to Paper I. 

The tree amplitude for beta decay of the polarised neutron is given by 

Up{P2)W X (p2,Pl)^ (l + S7 5 7 M ^) «n(Pl) , (1) 

where Gy = GfcosOq with the universal Fermi coupling Gf and the Gell-Mann-Levy- 
Cabibbo angle 6c, is the polarisation vector of the neutron with n 2 = —1, n ■ p\ = 0, and 
s = ±1, and W\(p2,Pi) is a general form of the weak vertex of hadrons and reads 

W x ( P2:Pl ) = 7 a(/v - 9A1 5 ) (2) 

1 It is included in Sirlin's proof [j] that the inner radiative corrections can be factored out for the processes 
concerning the unpolarised neutron. As far as the authors know there is no proof which shows that the inner 
radiative correction for the polarised neutron and that for the unpolarised neutron should agree. 
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at low energies. We retain fv = 1 to trace the vector coupling in the calculation. The 
spinors of the neutron, proton, electron and antineutrino are denoted by u n , u p , u e , and 
v u , respectively, with the momenta specified in parentheses. After spin summation and 
integration over p u the amplitude square reads 

£ \M^\ 2 = \QG 2 v m n m p E u [(/* + 3g\)E + 2(f v g A - g\)(n ■ £)] , (3) 

spin 

where E and E v are energies of the electron and the antineutrino. Therefore, the asymmetry 
parameter is given by 



A = " KJ ; 2 yA y 2 AJ , (4) 



2 (fygA - g\ ) 
fv + 3g\ 

the electron velocity factor j3 = \£\/E being removed as a convention. 

To evaluate full one loop corrections, we divide the integration region of the virtual gauge 
bosons into long- and short-distance parts 5 : 

(i) < \k\ 2 < M 2 , (ii) M 2 < \k\ 2 < oo, (5) 

where k is the momentum of the virtual gauge bosons, and the mass scale M, introduced by 
hand, divides the low- and high-energy regimes and is supposed to lie between the proton- 
neutron masses (m p and m n ) and the W and Z boson masses, m\y and mz- Old-fashioned 
four- Fermi interactions are applied to the proton and neutron in region (i), and the mass 
scale M is regarded as the ultraviolet cutoff of the QED (i.e., purely photonic) correction. 
In region (ii), electroweak theory is used for quarks and leptons, and M is the mass scale 
that describes the onset of the asymptotic behaviour. The concern is to connect the results 
in (i) and (ii) smoothly. Abers et al. |Hj proved that the logarithmic divergences that are 
proportional to f v are universal for the Fermi transition on the basis of the conserved vector 
current with the use of the current algebra technique. The same was proven for the g\ 
terms for the Gamow- Teller transition for which current conservation is broken only with 
soft operators [HI |H] ■ This guarantees smooth connection of the logarithmic divergence for 
the corrections of f v and g\. There appear, however, intereference terms of the order fvgA, 
for which logarithmic divergences depend on the model of hadrons. Marciano & Sirlin ^0] 
proposed the prescription to evaluate the high and low energy contributions separately by 
rendering the UV divergence in the low-energy contribution milder by taking account of form 
factors of hadrons. We follow the same prescription |1U1 II 1 j in the present calculation of the 
asymmetry parameter. 
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The diagrams of QED one- loop corrections are depicted in Figure 1, where (v) is the 
vertex correction, (s) is the self energy correction and (b) is bremsstrahlung. We write the 
one-loop amplitude 

M' = M {V) +M {s \ (6) 

The bremsstrahlung contribution is added separately. We consider the static limit for nucle- 
ons, q 2 = (pi — p 2 ) 2 <C m 2 . Our calculation is done in the Feynman gauge. 
We start with the vertex correction, which is given by 

d 4 k 1 11 



2^2 J (2vr) 4 (£ - k) 2 - m 2 e (p 2 + k) 2 -m 2 p k 2 - X 2 
xu e {l)r {7 • (£ - k) + m e } 7 A (1 - 7 5 )v u ( Pu ) 

x u p (p 2 hf 1 {7 • (P2 + k) + m p } W\(p 2 + fc,pi)(l + S7 5 7 • n)u n (p{), (7) 
where A is the photon mass to regulate the infrared divergence. Using identities, 

« e (*)7" {7 •(*-*:)+ »M = u e (£){(2t-ky + ia^K}, (8) 
u p {p2)l^ {7 • (P2 + k) + m p } = u p (p 2 ) {(2p 2 + k)fj, - ia pu k u } , (9) 

we decompose $7J) into three parts, 

M (v) =M (vi) +M W + M {v3) . (10) 

Here picks up the product of (2£ — kY in ((SJ) and (2^2 + A;)^ in ©, and at the same 

time W\(p2 + k,pi) is replaced by W\(j>2,Pi)- It has the same gamma matrix structure as 
the Born term If!]), and is then written as a multiplicative correction factor. This correction 
has both UV and IR divergences and depends on the electron velocity. The UV divergence 
in M^"" 1 ^ is cancelled by that in the self energy correction of Ai^ s \ 

= Uz 2 (m e ) - 1 + ^Z 2 (rn p ) - (11) 

and the IR divergence, along with that arising from MS 8 ', is cancelled when the contribution 
of bremsstrahlung M.^ is added. 

The term JH( v2 "> represents the combination of ia^ v k v in JSJ) and {2p2 + k)u in @, and 
W\(p 2 + k,pi) is again replaced with W\(p 2 ,pi). This term is UV and IR finite, but gives 
an electron-velocity dependent factor. In the static limit of the nucleon the correction is also 
a multiplication on the tree amplitude. 
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A straightforward calculation yields 



E {(-^ + A^ 2 W 0) * + cc.} + / e \M 

spin ^ ' spin 



(6) 



2 (Eq - E - oj) 
{Eq - E) 



lQGym n m p E v ^—2 



9(E) - I) (f v + 3g\) E 



2(g(E)--Uf v g A -gl)sl 



n 



(12) 



where g(E,Eo) is the conventional g function that appears in the radiative correction for the 
beta decay rate and is defined with an additional constant 3/4 |12l 17]. 



g(E,E ) 
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2(E Q - E) 



^-tanh- 1 /? I 2(1 + p A ) 4 



3£ 

2^ , (^o-i?) 2 



6£ 2 



4tanh" 1 ^l , 



(13) 



where Eq is the end point energy of the electron, and g(E,Eo) is a similar function for the 
spin-dependent term, 



g(E,E 



3 In 



in, 



in. 



3 + 4 l( 2/3 
A + p L [iTp 



1 



+4 ( -tann" 1 /? - 1 



E - E 3 (£ - £) 2 ln 2(£; - E) 



3£/3 2 



24£ 2 /? 2 



+^tanh _1 /3 (l - tanh" 1 /?) • 



Here 



l(z)= r^io g (i-t 

JO I 



(14) 



(15) 



is the Spence function. We also define g(E, Eq) with an additional constant 3/4 as a conven- 
tion. These are the outer corrections, which we may write 



Sout = ^9{E, Eq) 



jg(E,E ) 



(16) 



and agree with the formulae given by Shann and by Garcia and Maya [3] . They also agree 
with the expression derived by Yokoo et al. P] up to a constant. 
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We now consider the remaining term ^ 3 - ) , 

d 4 k 



2V2 J (2vr) 4 (£ - k) 2 - m 2 (p 2 + k) 2 - m 2 k 2 - A 2 
xu e (£) {(2£ - kY + ia^K) 7 A (1 - 7 5 H(P*) 

xu p (p 2 )R^x(P2,Pi, k)(l + S7 5 7.n)n„(pi), (17) 



where 



Rh\(Pi,P2, k) = (2p 2 + k)^ {W X {P2 + k,p 1 )-W x {p 2 ,Pi)}- icr^ u k u W x (p2 + k,p ± ) 
~ -ia llu k u Wx{p2 + k,pi) 

~ -ia^A^C/v - <m 7 5 ) , (18) 

in the approximation of the point nucleons. It is only this .M^ 3 ** term that depends on the 
details of the hadronic part of the weak current. It is clear from the powers of k that this 
term is IR convergent, whereas it is UV divergent. 
A straightforward calculation leads to 



spin 



= 16Gym n m p 



3 <?2 , , \ i M\ , /3, 2 . 9 



o (/v + fv9A) log — + 7 /v + 7 Ma > + * • 




3 ' 2 , , \ 1 / M\ 2 (1 2 , 5 



^(sA + Ma) log ^— J + + -Ma ) > (3EE U - ( ■ pi ) 



{(^U^g) 



g/y + gSU + ^Ma 



(E(n-p v ) + E v (n-£)) 



2 



+2* <j ^(ffi + Ma) log ( ^ 1 + r-g\ + \Sv9a) \ (E(n • pi) - E u (n • £ 



(19) 



We use the current algebra technique to classify logarithmic divergences into those that 
have universal coefficients irrespective of the model of hadrons and those that are model 
dependent. Repeating the same calculation as in Paper I but including the spin projection 
operator, we find that (3/2)/ y log(M/m p ) 2 in the first, (3/2)g\ \og(M/m p ) 2 in the second, 
(3/2) fygA log(M/m p ) 2 in the third and (3/2)g\ log(M/m p ) 2 in the fourth curly brackets are 
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universal. This observation tells us that the terms are summarised as, 



spin 



16G v m n m p 



1 
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4 8vr 2 
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-2sg 



2 I xGT' 

A ' in ' 4 8^ 2 



+ 



{^(n-^-^n-f)} 



(20) 



Adding the tree term and after integration over p v 

j2 \m^\ 2 + {m {v ^m {0 > + M {v3 >M 



(o)\ 



spin 



spin 



l6G v m n m p E u 1 + ^ • 



3 e" 

4 ' 8tt 2 



{/^ (l + CO + 3^ (l + 5g T/ ) } ^ 
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where 



5 F ' = — 
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3, (M r 
o lo S 



+ 



(21) 

(22) 
(23) 



are the inner corrections for Fermi and Gamow- Teller transitions that are defined in Paper 
I and the factor [1 + (3/4)e 2 /87r 2 ] is to be included in <5 out . In eqs. (J22J) and (J23J) the first 
logarithms are model- independent and the second with the coefficients gA/fv or fv/gA are 
model dependent. The correction from J\4^ v ^ is written as multiplicative factors on the 
coupling constants for both Fermi and Gamow- Teller parts while they are divergent within 
QED. 

The short distance correction from the integration region (ii) in (JSJ) is evaluated using 
electroweak theory 5 . When we consider corrections relative to muon decay, we only need 
to consider the box diagrams of photon (or Z) and W exchanges (see Fig. 2 of Paper I). In 
order to connect the quark-level amplitudes with hadronic ones, we assume that the ratio of 
the tree and loop amplitudes for beta decays of the assembly of quarks is the same as that 
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for neutron beta decay [3]. This is justified at least for the universal logarithmic divergent 
part. With this prescription the correction amounts to a multiplication factor 




3 5 
2 + 2tan 4 6»vy 




(24) 



where Q = 1/6 is the mean charge of the isodoublet of quarks, and the second line is the 
correction that appears in muon decay and thus subtracted when we consider the radiative 
correction relative to muon decay that determines Gf- 

This electroweak one-loop correction amounts to adding to 8f extra terms, 



Si = 




where the terms proportional to gA/fv are collected in C F , 




and similarly for , 



e 

8^2 



log 



m- 



+ l + 3Qlog 



m z 
M2 



+ c 



GT 



where C is 



(25) 



(26) 



(27) 



C 



GT 



9 A I 2 \ mi, 




(28) 



for point nucleons. 

We observe that the M dependence (upper cutoff) that appears in the first term of (|2*2*|) 
is cancelled by the first term in the braces in ([24)1 . which demonstrates a smooth connection 
from electroweak theory to effective hadronic theory for the Fermi transition. The UV di- 
vergence in the term proportional to gA/fv, however, fails to cancel against the divergence 
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with the coefficient 3Q, unless Q has a specific (and unrealistic) value of charge. Marciano 
and Sirlin ^Uj proposed to evaluate the model-dependent long-distance divergence of the 
Fermi transition by rendering it softer introducing nucleon form factors, and leave the term 
2>Q\og{raz / M) as it is taking M as the mass scale of the onset of the asymptotic behaviour 
|lUj . The same procedure was followed in Paper I for the GT part, while it was noted that the 
inclusion of the weak magnetism is important to evaluate the long-distance integral, especially 
for the GT transition. 

The calculation of the long distance contribution is made replacing 

rf ^ ^ Fl (fc2)_ * ^^(fc*), (29) 

2m n 

for the electromagnetic vertex, and 

W x (p 2 , Pl ) -» lx {f v F v (k 2 )-g A7 5 F A (k 2 )}--^-a^k u F w (k 2 ), (30) 

in eq. ([7jl. From the form we observed in eq. (|21j) we expect that the calculation incorporating 
form factors would give rise to a result summarised in the same form, while C F and C GT are 
modified exactly as in jSJ. Since we have not found an immediate proof that it should, we 
repeated a long calculation as we did in Paper I including the spin projection operator, and 
confirmed the anticipated result. In fact, we obtained C F and C GT exactly those that appear 
in the spin independent part. So we take the result of numerical integral of Paper I, 

C F = 1.751 + 0.409 = 2.160, (31) 
C GT = 0.727 + 2.554 = 3.281, (32) 

where the two parts of numbers represent contributions from the (V,A) interaction and weak 
magnetism. The first number in C F was evaluated by Marciano and Sirlin ^U] and by Towner 
jllj . and agrees with their results up to slight differences in the input parameters. 

In conclusion the radiative correction to polarised neutron beta decay to order 0(a) is 
summarised as 

d3 k ^..^(Eq-E-u) 



(2tt) 3 2u; ^ 1 1 (Eo-E) 

spin v ' spm v ' 

= l6G 2 v m n m p E u {/£ (l + 5f n + <U t ) + 9a (l + ^ T + *ant) } E 
+2s {f v g A (l + itf? + ^5g T + 4ut) - g\ (l + € T + <5out) } (ft ■ I) 



(33) 
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in the static nucleon approximation. Hence the asymmetry parameter is written as 

l + %g(E,E ) fvU 



A 



l + ^g(E,E ) f^ + 3g A 



(34) 



where 



ft 
A 



g 2 A (l + S\ 



in J 



(35) 



The denominator of eq. (|34j) is the combination that appears in the neutron decay rate. 
The energy dependent prefactor C(E), 



l + C(E) 



l + ^g(E,E ) 
Zn 



I l + —g(E,E Q ) (36) 

is plotted in Figure 2 as a function of the kinetic energy T = E — m e . The magnitude of 
(a/2n)g(E, Eq) and (a/2ir)g(E, Eq) is about 2%, but the two corrections nearly cancel in 
g(E, Eq) — g(E, Eq), leaving the net outer correction for the asymmetry being quite small, of 
the order of 0.1%. For convenience we give a fit to C(E) for neutron beta decay with 



C(E) = -0.00163 + 0.00210/£ + 0.000491£, 



(37) 



where E is in units of MeV. The fit, also displayed in Figure 2, overlays nearly top on the 
true function of C(E). 

The cancellation also takes place for the inner correction. After correcting for C(E), the 
axial-vector to vector coupling ratio extracted from the tree level formula is related to its 
tree-level value as 



= 1.0012 (38) 
JV 

The dominant part of the inner correction, including log mz/m p cancels in 5f n — and the 
net correction is of the order of 0.1% for gA/fv (which is usually denoted as A = Ga/Gv = 
-QA/fv)- 

The Particle Data Group ^3] gives a value gA = 1.2670 ± 0.0030. This is obtained by 
averaging 5 values reported in the literature, one of which ^1] is obtained including the 
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outer radiative correction (with the inner correction discarded), and others are results that 
do not include radiative corrections. The outer radiative correction reduces the value of 
\gA\ by about 0.0007, but the scatter among the data from different authors is 0.005 (rms), 
so systematic errors other than the radiative correction dominate the uncertainty of gj±. 
As we have shown that the inner radiative corrections can be included into qa i n common 
irrespective of quantities measured for beta decay, it is a matter of definition whether they 
are included in g& or not. If we define the tree-level axial coupling constant it is related with 
the value including the radiative correction by eq. ((3*5|) . 
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Figure Captions 



Fig. 1 radiative corrections to neutron beta decay. 

Fig. 2 Outer radiative correction C(T + m e ) for the asymmetry parameter as a function of 
the kinetic energy of electron. The solid curve is C(T + m e ), and the dotted curve, which 
overlays nearly exactly on the solid curve, is fit (37). 
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